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In addition to being critical for apoptosis, components 
of the apoptotic pathway, such as caspases, are in- 
volved in other physiological processes in many types 
of cells, including neurons. However, very little is known 
about their role in dynamic, nonphysically destructive pro- 
cesses, such as axonal arborization and synaptogenesis. 
We show that caspases were locally active in vivo at the 
branch points of young, dynamic retinal ganglion cell 
axonal arbors but not in the cell body or in stable mature 
arbors. Caspase activation, dependent on Caspase-3, 
Caspase-9, and p38 mitogen -activated protein kinase 



(MARK), rapidly increased at branch points correspond- 
ing with branch tip addition. Time-lapse imaging revealed 
that knockdown of Caspase-3 and Caspase-9 led to more 
stable arbors and presynaptic sites. Genetic analysis 
showed that Caspase-3, Caspase-9, and p38 MARK in- 
teracted with Slitla-Robo2 signaling, suggesting that lo- 
calized activation of caspases lie downstream of a ligand 
receptor system, acting as key promoters of axonal branch 
tip and synaptic dynamics to restrict arbor growth in vivo 
in the central nervous system. 



Introduction 

The formation of functional synaptic connections is critical for 
the normal operation of the nervous system. In the central ner- 
vous system (CNS), synapse formation depends on the regulated 
elaboration of the arbors (consisting of branches) of pre- and 
postsynaptic cell axons and dendrites, processes that can occur 
concurrently (Waites et al., 2005; Meyer and Smith, 2006; 
Ruthazer et al., 2006; Jin and Gamer, 2008; Gibson and Ma, 
2011). The processes of arborization and synaptogenesis are cru- 
cial in defining neuronal morphology and synaptic connectivity 
(Lin and Koleske, 2010). Abnormalities in arbor structure can 
lead to changes in behavior (Smear et al., 2007) and neurodevel- 
opmental disorders (Wong, 2008). Arbor and synaptic structure 
are often affected early in the progression of neurodegenerative 
conditions (Wishart et al., 2006; Saxena and Caroni, 2007), il- 
lustrating the fundamental importance of the regulation of these 
processes in the nervous system. 

The vertebrate retinotectal projection has become a power- 
ful model system for the study of axon pathfinding, arborization, 
and synaptogenesis (Dingwell et al., 2000; Hutson et al., 2004; 
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Cohen-Cory, 2007; Poulain et al., 2010). The growth cones at 
the tips of retinal ganglion cell (RGC) axons navigate along a 
stereotypical pathway from the retina to their target, the optic 
tectum in the midbrain. Upon reaching the tectum growth cones 
change their morphology with the production of branches to 
develop an arbor, accompanied by the formation of synapses, 
both of which are highly dynamic processes (Harris et al., 1987; 
Meyer and Smith, 2006; Ruthazer et al., 2006; Campbell et al., 
2007). A variety of ligand receptor systems have been identified 
taking part in these processes, including brain-derived neuro- 
trophic factor (Cohen-Cory and Fraser, 1995; Alsina et al., 
2001), Netrin-1 (Manitt et al., 2009), SemaSA (Semaphorin 3A; 
Campbell et al., 2001), and Slitla-Robo2 (Roundabout 2) sig- 
naling (Campbell et al., 2007). However, unlike the intensive 
investigation of the intracellular signaling pathways acting 
downstream of ligand receptor systems operating in growth cone 
guidance (O'Donnell et al., 2009; Bashaw and Klein, 2010), our 
understanding of the intracellular signaling pathways regulating 
arborization and synaptogenesis remains much more limited. 

Neurons are compartmentalized into a cell body, an axon, 
and dendrites. Somaless Xenopus laevis axonal growth cones 
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Figure 1 . Local caspase activation at branch points in young RGC arbors. (A, D, E, and H-R) Representative IMD ratio images of 80-hpf live RGC arbors 
in the optic tectum (A, D, E, and H-P) and an RGC cell body (Q and R) in the retina. Caspase activation is represented by the pseudocolors that correspond 
to the Venus/ECFP ratio (1 .5-0.5). Red represents low caspase activation, and violet represents high caspase activation. (A, D, E, and H-P) Control 
(A and D), SCAT3 (E and H), SCAT3 + coMO (I and J), SCAT3 + cosp-3aM01 (K and L), SCATS + casp-9M01 (M and N), and SCAT3 + p38mapkMO] 
(O and P) are shown. (B-D and G) ECFP channel (B and D) and Venus channel (C and G) of SCAT3- and control-expressing arbors ore also shown. 
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are able to undergo chemotropic responses in vitro (Campbell 
and Holt, 2001; Ming et al., 2002) and to navigate to and arbo- 
rize in the optic tectum in vivo (Harris et al., 1987), suggesting a 
role for local signaling pathways. Indeed, the structural com- 
partmentalization of neurons allows the partitioning of cellular 
signaling pathways, such as those regulating protein synthe- 
sis and degradation to axons and/or dendrites (Steward and 
Schuman, 2003; Cajigas et al, 2010; Bingol and Sheng, 2011; 
Jung et al., 201 1). This could allow very precise local control of 
neuronal morphology and remodeling critical for the correct de- 
velopment and function of neurons. Caspases are an example of 
a protease system, originally identified for their critical roles in 
apoptosis and cytokine processing (Creagh and Martin, 2001; 
Ulukaya et al., 201 1). More recently, localized nonapoptotic roles 
for caspases have been identified (Kuranaga and Miura, 2007; 
D'Amelio et al, 2012; Hyman and Yuan, 2012). For example, in 
neurons, caspases may be active and act in axons or dendrites, 
mediating the chemotropic responses of Xenopus retinal growth 
cones (Campbell and Holt, 2003), dendritic pruning during 
metamorphosis in Drosophila melanogaster (Kuo et al., 2006; 
Williams et al., 2006), mammalian long-term depression (Li 
et al., 2010; Jiao and Li, 201 1), and songbird learning (Huesmann 
and Clayton, 2006). However, the role of the apoptotic pathway 
and the principles underlying its regulation in local nonapoptotic 
roles in neurons in vivo, such as during the processes of axonal 
arborization and synaptogenesis, remain largely unknown. 

To address these issues, we imaged caspase activity in the 
CNS in vivo, in the intact living zebrafish embryo. We show 
that caspase activity is very local, at the branch points of young 
dynamic RGC axonal arbors, in a Caspase-3-, Caspase-9-, and 
p38 MAPK-dependent manner, but not in RGC cell bodies. 
Time-lapse imaging revealed that caspase activity undergoes 
rapid activation at branch points during branch tip formation. 
In contrast, older, more stable arbors have reduced levels of 
caspase activation. Furthermore, a reduction of Caspase-3, 
Caspase-9, and p38 MAPK function leads to more stable arbors 
and interactions with Slitla-Robo2 signaling to restrict arbor 
growth and synaptogenesis. Our findings suggest that local cas- 
pase activation is a key promoter of axon remodeling in vivo in 
the vertebrate CNS. 

Results 

Local caspase activafeion at branch points 
of young RC3C arbors 

To investigate whether caspases play nonapoptotic roles in the 
processes of arborization and synaptogenesis, we performed 
mRNA in situ hybridization at 80 h postfertilization (hpf) with 



zebrafish casp-3a (caspase-Sa) and casp-9 {caspase-9), key 
effector and initiator caspases of the intrinsic apoptotic pathway, 
and casp-8 (caspase-8), an activator of Caspase-3 in the extrinsic 
apoptotic pathway, respectively (Ulukaya et al., 2011). casp-Sa 
and casp-9 are strongly expressed in the RGC layer of the retina 
and throughout the tectum (Fig. SI, A-D), whereas casp-8 
showed a similar though much weaker expression pattern (Fig. SI, 
E and F). Next, we used the caspase activity reporter, SCAT3 
(sensor for activated caspases based on fluorescence resonance 
energy transfer [FRET] 3; Takemoto et al., 2003, 2007) to reveal 
the temporal and spatial pattern of caspase enzymatic activity in 
the retinotectal system in vivo. SCAT3 consists of ECFP and 
Venus separated by a peptide linker region containing a DEVD 
(Asp-Glu-Val-Asp) consensus Caspase-3 cleavage site. The 
DEVD cleavage site specificity is conserved for zebrafish 
Caspase-3 (Yabu et al., 2001). Because the peptide linker region 
of SCAT3 is cleaved by caspase activity, its activity is reported 
as a decrease in the Venus/ECFP ratio, representing a decrease 
in FRET (Takemoto et al., 2003). SCAT3 or control constructs 
lacking the cleavage site were transiently expressed in RGCs 
and tectal cells (Materials and methods; Campbell et al., 2007; 
Sato et al., 2007). SCAT3-expressing 80-hpf RGC axonal ar- 
bors showed a significant decrease in the Venus/ECFP over the 
whole arbor and at branch points compared with surrounding 
regions (regions of interest) defined in the legend to Fig. 1 , rela- 
tive to arbors expressing the control construct (Fig. 1, A-H and 
quantified in S-U). 

To demonstrate that the decrease in Venus/ECFP ratio 
observed for SCAT3-expressing arbors primarily reflected 
Caspase-3 activity, we used antisense morpholino (MO) oli- 
gonucleotides to partially knockdown Caspase-3 function. A 
previously published MO against zebrafish casp-3a {casp- 
JaMOl; Yamashita et al., 2008) blocked the decrease in the 
SCAT3 Venus/ECFP ratio over the whole arbor and at branch 
points compared with surrounding regions, whereas a five-base 
pair mismatch control MO (controlMO [coMO]) did not (Fig. 1, 
I-L and quantified in S-U). Additionally, knockdown of casp-8 
(casp-SMOl; Sidi et al., 2008), whose activity may partially 
cleave SCAT3 in vitro (Takemoto et al., 2003), did not (Fig. 1, 
T and U). Because it has not previously been reported that 
SCAT3 is an effective reporter of Caspase-3 activity in the 
zebrafish, we performed additional control experiments (Fig. S2). 

To further confirm apoptotic pathway activation of SCAT3, 
we knocked down casp-9, a key initiator caspase, which proteolyti- 
cally cleaves and activates executioner caspases, such as Caspase-3 
(Ulukaya et al., 2011), and whose activity may partially cleave 
SCAT3 in vitro (Takemoto et al., 2003). Compared with coMO 
morphant arbors, knockdown of Caspase-9 function with two 



The region enclosed by the dashed squares in A, E, I, K, M, O, and Q ore magnified in B-D, F-H, J, L, N, P, and R. (S) Quantification of Venus/ECFP ratios 
across the whole arbor. Branch points were determined as covering a distance of 1 pm or less away from where a branch tip joins an arbor. Surrounding 
regions were determined at a distance of 1 pm away from a branch tip and 2 pm in length, covering the width of the branches. (B-D and F-H] Examples 
illustrating regions of interest for the quantification of the Venus/ECFP ratio at branch points (purple dashed areas) and surrounding regions (white dashed 
areas) for control- and SCAT3-expressing arbors, respectively. (T) Quantification of Venus/ECFP ratios at branch points (BP) and surrounding regions (S). 
(U and V) Presented as a ratio of branch points/surrounding regions (U) and RGC cell bodies (V). White asterisks indicate the parent axon. 19-28 RGC 
arbors were analyzed in S. 29-36 branch points and 87-1 08 surrounding regions were analyzed in T and U. 1 5-1 6 RGC cell bodies were analyzed in V. 
Error bars represent SEMs. *, P < 0.05; **, P < 0.01 ; ***, P < 0.001 . Bars, 1 0 pm. The image in E of a SCAT3-expressing arbor is shown again in Figs. 5 A 
and S2 C. The magnified image in H also is shown again in Fig. 5 B. 
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Figure 2. Caspase activity is undetectable in older RGC arbors. |A-D) Representative IMD ratio images of 100-hpf live RGC arbors in the optic tectum; 
caspase activation is represented by the pseudocolors that correspond to the Venus/ECFP ratio (1 .5-0.5). Red represents low caspase activation and violet 
represents high caspase activation. (A-D) Control (A and B) and SCAT3 (C and D) images are presented. The region enclosed by the dashed squares in 
A and C are magnified in B and D. (E-G) Quantification of Venus/ECFP ratios across the whole arbor for both 80-hpf (E; data as in Fig. 1 S) and 1 00-hpf 
(F) arbors at branch points (BP) and surrounding regions (S) presented as a ratio of branch points/surrounding regions (G). White asterisks indicate the 
parent axon. Dorsal views are shown, and anterior is up. 10-12 RGC arbors were analyzed per condition. Error bars represent SEMs. **, P < 0.01 . 
Bars, 1 0 pm. 



different MOs (translation-blocking casp-9M0l and splice- 
blocking [SB] casp-9SBM0l ; Sidi et al., 2008) blocked the de- 
crease in Venus/ECFP ratio (Fig. 1, M and N and quantified in 
S-U), indicating that SCAT3 reports caspase activation. Cas- 
pase activation was also blocked by knockdown of p38 MAPK, 
expressed in zebrafish RGCs (Krens et al., 2006), and which we 
have shown previously to be upstream of Caspase-3 activation 
in retinal growth cones (Fig. 1, O and P and quantified in T and U; 
Campbell and Hok, 2003; Hsu et al., 201 1 \p38mapkMOi]). As 
there is a low amount of apoptotic cell death in zebrafish RGCs 
at 80 hpf (Biehlmaier et al., 2001), we next tested that the lack 
of apoptosis observed is caused by lack of caspase activation in 
RGC cell bodies. SCAT3- or control-expressing 80-hpf RGC 
cell bodies did not show caspase activation (Fig. 1, Q and R 
and quantified in V). We were also unable to detect significant 
caspase activation in the dendrites of 80-hpf and 100-hpf tectal 
cells (Fig. S3). 

Caspase activity rapidly increases upon 
branch tip formation 

Caspase activation at branch points in 80-hpf RGC arbors corre- 
lates with a period during which they are highly dynamic 
(Stuermer, 1988; Meyer and Smith, 2006; Campbell et al., 
2007). Older (100 hpf) arbors that we have previously shown to 
be more stable (Stuermer, 1988; Campbell et al., 2007) had a 



high Venus/ECFP ratio, i.e., low caspase activation (Fig. 2, A-D 
and quantified in E-H). To establish whether caspase activation 
is dynamically regulated, we visualized branch point additions 
and retractions directly over time in the same arbors, via time- 
lapse imaging of 80-hpf and 100-hpf arbors expressing SCAT3 
every 5 min for ^2 h. For 80-hpf arbors, we observed branch tip 
additions and measured the Venus/ECFP ratio for the region of 
the arbor that would become a branch point in the next 5 min and 
a neighboring region of the same arbor that would not, i.e., dy- 
namic versus nondynamic areas of the same arbor (Fig. 3 A, re- 
gion of interest defined in the legend and imaging scheme). We 
found that a region of an arbor that became a branch point in 
5 min underwent a significant increase in caspase activation 
during the 5-min interval between time points compared with a 
neighboring region of an arbor that did not become a branch 
point (Fig. 3, B-E and quantified in F and G; and Video 1). Im- 
aging of stable branch points and neighboring regions of 100-hpf 
arbors, which do not have a branch point, nor will they develop 
one during the 5-min imaging interval (Fig. 3 H), did not show 
significant caspase activation (Fig. 3, I-L and quantified in M 
and N; and Video 2). Thus, caspase activity rapidly increases 
(within 5 min) at branch points, immediately before or during 
branch tip formation, whereas a reduction of caspase activation 
directly correlates with increased branch tip stability. We were 
unable to detect significant caspase activation at branch points 
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Figure 3. Caspase acHvity rapidly increases 
upon young RGC arbor branch tip formation. 

(A) Imaging analysis scheme for 80-hpf ar- 
bors expressing SCAT3-forming branches dur- 
ing the 5-min imaging interval. Branch points 
were determined as covering a distance 
of 2 ^jm or less away from where a branch 
tip joins, and arbor nonbranch points were 
5-20 pm away from a branch tip and 3 (jm 
in length, covering the width of the branches 
as for surrounding regions (Fig. 1). (B-E) Rep- 
resentative IMD images of an 80-hpf RGC 
arbor expressing SCAT3. Caspase activa- 
tion is represented by the pseudocolors that 
correspond to the Venus/ECFP ratio 1 .5-0.5 
or 1.0-0.5. Red represents low caspase ac- 
tivation, and violet represents high caspase 
activation. (B-D) 5 min before branch point 
formation (BP T — 5 min; B) and after branch 
tip [BP] formation (C). The region enclosed by 
the dashed squares in B and C are magni- 
fied in D and E. (F and G] Quantification of 
the Venus/ECFP ratio (F) comparing changes 
in arbors between regions that do or do not 
lead to the formation of branch points as in 
the scheme in A and presented as ratios (G). 
(hi) Imaging analysis scheme for 100-hpf ar- 
bors expressing SCAT3 that are stable during 
the 5-min imaging interval. (I-L) Representative 
IMD images of a 1 00-hpf RGC arbor express- 
ing SCAT3 that possesses branch points (BP 
TO; I), which remain stable during the 5-min 
imaging interval (BP T 5 min; J). The regions 
enclosed by the dashed squares in I and J are 
magnified in K and L. (M and N) Quantifica- 
tion of the Venus/ECFP ratio (M] comparing 
changes in arbors between stable branch tips 
and regions that do not form branch tips as 
in the scheme in H and presented as ratios 
(N). White arrowheads depict the region in D 
that will become a branch point in E. White 
asterisks indicate the parent axon. 37 branch 
point additions and nonbranch points in six 
arbors were analyzed for 80-hpf arbors. 30 
stable branch points from seven arbors were 
analyzed for 100-hpf arbors. Dorsal views 
are shown, and anterior is up. Error bars rep- 
resent SEMs. ***, P < 0.001 . Bars, 10 pm. 



where branch tips had undergone recent retractions (i.e., within 
5 min; Fig. S4). However, we were able to detect a significant 
decrease in caspase activation between regions of an arbor that 
underwent recent branch point retractions compared with those 
that had undergone recent branch tip additions (Fig. S4 H). 

Caspase-3, Caspase-9, and p38 MAPK 
promo'te RGC axonal arbor and presynaptic 
dynamics in vivo 

Our in vivo imaging studies using SCAT3 indicate that caspase 
activation is more active before and/or dming branch tip formation. 



suggesting that it may promote arbor dynamics. We performed 
time-lapse imaging every 15 min for 5 h of 80-hpf arbors ex- 
pressing the presynaptic marker YFP-Rab3 (Campbell et al., 
2007), taking advantage of our observations that casp-3aMOl, 
casp-9MOl, and p38mapkMOl were able to reduce caspase ac- 
tivation (Fig. 1, G-J and quantified in M-0), that there is no in- 
crease in apoptosis from the administration of the MOs, and that 
the RGC axons of the morphants project normally to the optic 
tectum (Fig. S5). casp-3aM0l, casp-9M0l, and p38mapkMOl 
morphant arbors showed a reduction in dynamic branch tips 
accompanied by an increase in stable branch tips, decreased 
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Figure 4. Capase-3, Caspase-9, and p38 
MAPK promote RGC axon arbor dynamics. 

(A, C, E, G, and I) Tracings of arbors at 80-85 
hpf expressing YFP-Rab3. Images represent 
WT (A), coMO (C), casp-3aMO] (E), casp- 
9M01 (G), and p38mapkl\AO] (I) morphant 
arbors every hour during a 5-h time-lapse im- 
aging. All show dorsal views. Block, stable 
branch tips; red, eliminated tips; green, added 
tips; blue, tips that were added then eliminated 
during imaging. (B, D, F, H, and J). The pro- 
portion of branch tips in each category. These 
graphs ore derived from analysis of a single 
arbor in a single experiment. (K and L] Branch 
tip additions (K) and retractions (L] ore shown 
in the graphs. (N and O) Presynaptic puncto 
gain (N) and loss (O) ore presented as a per- 
centage of initial branch tip or presynaptic 
puncto, respectively. (M and P] Total branch 
tip (M) and presynaptic puncta lifetimes (P) 
include all branch tips or presynaptic puncta 
present at one or more time points. (Q and R) 
Total branch tip (Q] and presynaptic puncta 
lifetime distribution (R] are shown. 9-1 3 RGC 
arbors were analyzed per condition. Tracings 
ore dorsal views; anterior is up. Error bars 
represent SEMs. *, P < 0.05; **, P < 0.01 ; 
***, P < 0.001. Bar, 10 pm. 
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branch tip additions, and retractions to approximately one 
third of those of wild-type (WT) or coMO arbors (Fig. 4, A-L; 
and Videos 3, 4, 5, 6, and 7). This was accompanied by an ap- 
proximate doubling in total branch tip lifetimes, a significant 
decrease in short-lived branch tips (with a lifetime of <15 min), 
and an increase in the proportion of branch tips present 



throughout the whole imaging period (Fig. 4, M and Q; and 
Videos 3, 4, 5, 6, and 7). Similar effects were observed for 
presynaptic puncta gain, loss, total puncta lifetimes, the pro- 
portion of short-lived puncta, and the proportion of puncta 
present during the entire imaging period (Fig. 4, N-P and 
R; and Videos 3, 4, 5, 6, and 7). Our findings indicate that 
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Figure 5. Caspase activation is reduced in RGC arbors with impaired Slit1a-Robo2 signaling. (A-F) Representative IMD ratio images of 80-hpf live RGC 
arbors in the optic tectum expressing SCAT3. Caspase activation is represented by the pseudocolors that correspond to the Venus/ECFP ratio (1 .5-0.5). 
Red represents low caspase activation, and violet represents high caspase activation. WT (A and B), ast (C and D), and slitla transMOl (E and F) are 
shown. The regions enclosed by the dashed squares in A, C, and E are magnified in B, D, and F. (G and H) Quantification of Venus/ECFP ratios (G) at 
branch points (BP) and surrounding regions (S) presented as a ratio of branch points/surrounding regions (H). White asterisks indicate the parent axon. 
15-35 branch points and 45-105 surrounding areas were analyzed per condition. Dorsal views are shown, and anterior is up. Error bars represent 
SEMs. *, P<0.05; **, P< 0.01 ; ** *, P < 0.001 . Bars, 1 0 pm. The image in A of a SCATS-expressing arbor is shown again in Figs. 1 E and S2 C. The 
magnified image in B also is shown again in Fig. 1 H. 



Caspase-3, Caspase-9, and p38 MAPK are key promoters of 
RGC axon arbor dynamics. 

Caspa5e-3, Caspase-S, and p3B MAPK 
interact with 5lit:1 a-RoboS signaling to 
restrict arbor growth 

The knockdown of Caspase-3 and Caspase-9 function leads 
to more stable branch tips and presynaptic sites, which is very 
similar to loss of Slitla-Robo2 signaling on RGC arbor dynam- 
ics (Campbell et al., 2007). This observation in combination 
with the observations that Robo2 acts cell autonomously in 
RGCs (Campbell et al., 2007) and that the apoptotic pathway is 
locally activated in RGC arbors leads to a possible model in 
which Caspase-3 and Caspase-9 might be part of the mecha- 
nism downstream of Slitla-Robo2 signaling, promoting arbor 
dynamics and restricting growth. To examine this, we expressed 
SCAT3 and observed caspase activity in 80-hpf astray {ast)/ 
robo2 mutant, slitla translation-blocking (trans) MOl {slitla 
transMOl), and morphant arbors at a dose known to increase 
arbor size and complexity (Campbell et al., 2007). Compared 
with age-matched WT or control morphant arbors, we were un- 
able to detect an increase in caspase activation at branch points 
relative to surrounding regions (Fig. 5, A-F and quantified in 
G and H), proposing that a reduction of Slitla or Robo2 



function leads to the earlier down-regulation of caspase activa- 
tion in RGC arbors and that caspase activation may act in the 
Slitla-Robo2 signaling pathway. 

As a functional test of the aforementioned model, we ex- 
amined gene-dosage effect interactions between astlrobol, 
slitla, casp-3a, and casp-9. We observed 80-hpf YFP-Rab3- 
expressing RGC arbors from WT, ast, astl+, robol, or slitla 
morphant embryos injected with either nothing, coMO, casp- 
JaMOl, or casp-9y[0\. Loss of one copy of Robo2 in ast/+ or 
knockdown of robo2, slitla, casp-Sa, or casp-9 led to arbors 
that were indistinguishable to WT or coMO arbors by all our 
quantitative parameters {astl+, robolSWsAO, slitla transMOl 
[low dose], cfljp-JflMOl, and casp-9MO\ compared with WT 
and coMO; Fig. 6, A-G and quantified in M-R). However, 
combining the reduction of Robo2 or Slitla function with ei- 
ther knockdown of casp-3a or casp-9 led to significant increases 
in branch tip number, arbor area, total arbor length, arbor 
order, and presynaptic puncta number, which were similar in 
size and complexity to ast mutant arbors {astl+; casp-3aMO\, 
robo2SBMO; casp-9M0\, slitla transMOl [low dose]; casp- 
3aM01, slitla transMOl [low dose]; casp-9M0\, and ast 
compared with control conditions; Fig. 6, H-L and quantified in 
M-R; Campbell et al., 2007). Combining the knockdown of 
p38mapk with the knockdown of slitla or robol also led to an 
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increase in arbor size, complexity, and presynaptic puncta number 
(Fig. 7, A-F and quantified in G-L). These observations implicate 
the apoptotic pathway and p38 MAPK as negative regulators of 
arbor growth via an interaction with SUtla-Robo2 signaling. 



To test for a genetic interaction between slitla, robo2, and 
casp-Sa in the regulation of branch tip dynamics, we performed 
time-lapse imaging for 5 h of 80-hpf YFP-Rab3-expressing 
RGC arbors using combinations of doses of MOs (low dose), 
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which on their own did not affect arbor dynamics (Fig. 8, A-F 
and K-M; and Videos 8 and 9). Combining robo2SBM0 or 
slit la transMOl with casp-SaMOl led to a reduction in dynamic 
and an increase in stable branch tips (Fig. 8, G-J; and Video 10). 
Significant decreases in branch tip additions and retractions were 
observed (Fig. 8, K and L), which were accompanied by a signifi- 
cant increase in mean branch tip lifetimes and the proportion of 
stable branch tips during the 5-h imaging period (Fig. 8, M and N). 
Furthermore, presynaptic puncta gain and loss were significantly 
reduced, accompanied by an increase in puncta lifetimes (Fig. 8, 
O-R). Collectively, these experiments imply a close genetic inter- 
action between Slitla-Robo2 signaling and the apoptotic pathway 
in the restriction of arbor, presynaptic growth, and in the promotion 
of arbor dynamics. 

□iscussion 

A role for the apoptotic pathway in 
promoting axonal arbor dynannics 

We have used a combination of in vivo imaging and genetics and 
found that the apoptotic pathway promotes axonal arborization 
and presynaptic dynamics in vivo, which are based on remodel- 
ing in terms of branch tip addition and retraction, rather than 
being physically destructive or degenerative. We identify novel 
roles for the apoptotic pathway in a local nonapoptotic context, a 
recently identified cellular phenomenon that is poorly under- 
stood. Arbor dynamics occur in response to attractive and repul- 
sive cues enabling pre- and postsynaptic sites to find appropriate 
partners with which to form stable synapses (Jin and Garner, 
2008). The promotion of axonal arbor and presynaptic puncta 
dynamics indicates that the apoptotic pathway in RGC axonal ar- 
bors promotes characteristics of young arbors, rather than matu- 
ration as in olfactory neurons (Ohsawa et al., 2010). 

Time course of caspase activation 

Our study is the first to examine caspase activation in a non- 
apoptotic role in vivo in the intact living organism. We found 
that caspase activation occurs rapidly within a 5-min period co- 
inciding with RGC branch point addition in vivo. This is consis- 
tent with our earlier in vitro data in which chemotropic cues 
were able to induce pro-Caspase-3 cleavage and enzymatic ac- 
tivation in 5 min (Campbell and Holt, 2003) and in which 
cleaved Caspase-3 was detectable after 10 min during the in- 
duction of long-term depression (Li et al., 2010). Because 
SCAT3 can be partially cleaved in vitro by Caspases-8 and -9 
(Takemoto et al., 2003), a small proportion of the decrease in 
FRET observed in RGCs may not be caused by Caspase-3 acti- 
vation. As the FRET ratio in casp-Sa morphant RGC arbors is 
indistinguishable to WT arbors expressing the control uncleav- 
able construct and knockdown of casp-8 does not affect the 
FRET ratio (Fig. 1), this possibility is less likely. As SCAT3 is 
cleaved upon caspase activation and the axonal cytoskeleton 
may restrict the diffusion of cytoplasmic proteins (Kulkarni 
et al., 2007), we were unable to fully assess the time course for 
the activation. Restiiction of the diffusion of proteins, such as 
caspases, in axons could be a mechanism for maintaining their 
activation locally. 



Relationship between caspase activation 
and upstream regulators 

Localization of caspase activity away from the cell body may 
represent an underlying principle of nonlethal apoptotic path- 
way activation in neurons, complementing those of lower 
sublethal and transient caspase activation (Jiao and Li, 2011). 
The upstream pathways regulating local caspase activation in 
neurons are largely unknown. We show here that the caspase 
activation reported by SCAT3 is dependent on Caspase-3, 
Caspase-9, and p38 MAPK. Because the effect of knockdown 
of Caspase-3, Caspase-9, and p38 MAPK function on RGC dy- 
namics is similar to the decrease in branch tip and presynaptic 
puncta dynamics observed in the loss of Slitla-Robo2 signal- 
ing (Campbell et al., 2007), we tested for and identified an 
interaction between Slitla-Robo2 signaling, p38 MAPK, and 
caspase activation and function. In our previous findings that 
Robo2 acts cell autonomously in RGCs, slitla is strongly ex- 
pressed in the tectum (Campbell et al., 2007), and here, the fact 
that caspase activation occurs at the branch points of RGCs 
suggests a potential model in which Slit la-induced Robo2 sig- 
naling in RGCs via p38 MAPK leads to local apoptotic path- 
way activation, possibly via localized retention of Slit protein 
by the extracellular matrix of the tectum (Xiao et al., 2011). 
The promotion of arbor dynamics and restraining of arbor 
growth is consistent with preventing the premature maturation 
of RGC arbors. Other attractive candidates for ligand receptor 
systems upstream of caspase activation include known acti- 
vators, such as L-a-lysophosphatidic acid (Steiner et al., 2000; 
Campbell and Holt, 2003), Sema3A (Shirvan et al., 1999), and 
Netrin-1 (Campbell and Holt, 2003). Netrin-1 is particularly 
interesting because it is able to stimulate Caspase-3 activation 
in vitro (Campbell and Holt, 2003) and promote RGC arbor 
dynamics in vivo (Manitt et al., 2009), though the relationship 
between the two remains unexplored. 

Knockdown of casp-9 leading to increased arbor stability 
is consistent with a role for the intrinsic apoptotic pathway 
(Ulukaya et al., 201 1). This finding is in contrast to our earlier 
in vitro study (Campbell and Holt, 2003) in which we were un- 
able to identify a role for Caspase-9 in chemotropic turning. A 
potential reason could be a result of the choice of guidance 
cues used (Netrin-1 and lysophosphatidic acid; Campbell and 
Holt, 2003), compared with the interactions of multiple ligand 
receptor pathways, which could operate in vivo to regulate 
local caspase activation. The role of other components of the 
intrinsic pathway, such as mitochondria, present in axons un- 
dergoing extensive dynamics (Zinsmaier et al., 2009), remains 
to be determined. A small proportion of SCAT3 cleavage may be 
caused by Caspase-9 activity (discussed in the first paragraph 
of the Discussion), and it remains possible that Caspase-9 also 
functions independently of Caspase-3. Additionally, Caspases-6 
and -8 are present in axons (Carson et al., 2005; Nikolaev et al., 
2009), and they may be potentially involved in arbor remod- 
eling independently of Caspases-3 and -9. Furthermore, be- 
cause very little is known about the presence, localization, 
and function of the broader caspase family in axons, it re- 
mains an open question whether other members may func- 
tion in arbor remodeling. 
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Figure 8. Caspase-3 interacts with Slit1a-Robo2 signaling to promote arbor dynamics. (A, C, E, G, and I) Tracings of arbors at 80-85 hpf expressing 
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Possible links to the cytoskelefeon and 
axonal degeneration 

Axonal branch tip and presynaptic puncta dynamics are ulti- 
mately regulated by the underlying cytoskeleton (Jin and 
Garner, 2008; Gallo, 2011). Caspase activity regulates cyto- 
skeletal dynamics during apoptotic (Ndozangue-Touriguine 
et al., 2008) and nonapoptotic processes (Kuranaga et al., 2006; 
Oshima et al., 2006) and cleaves numerous substrates (Fischer 
et al., 2003; Timmer and Salvesen, 2007; Crawford and Wells, 

2011) . Substrates important for axonal arbor remodeling may 
include actin, actin-binding proteins, and those in signal trans- 
duction pathways regulating cytoskeletal dynamics, for exam- 
ple, gelsolin, an actin-severing and -capping protein (Kothakota 
et al., 1997). Interestingly, cultured hippocampal neurons from 
gelsolin knockout mice show decreased filopodial dynamics 
(Lu et al., 1997), reminiscent of the effect of loss of Caspase-3 
function on arbor dynamics. 

Local apoptotic pathway activation via caspases plays 
roles in the destructive processes of developmental pruning and 
axonal degeneration (Kuo et al., 2006; Williams et al., 2006; 
Nikolaev et al., 2009; Schoenmann et al., 2010; Simon et al., 

2012) . Our findings that the apoptotic pathway promotes axonal 
arbor branch tip dynamics may indicate that an optimum amount 
of or localization of caspase activity is necessary for these non- 
destructive remodeling processes to occur. Because changes in 
synaptic function, arbor structure, and axonal degeneration 
occur early during neurodegenerative conditions, such as Al- 
zheimer's and Parkinson's diseases (Wishart et al., 2006; 
Saxena and Caroni, 2007; D' Amelio et al., 201 1, 2012), identi- 
fying the mechanisms of local apoptotic pathway activation in 
axons may be relevant to understanding these disorders. We 
propose that the apoptotic pathway is critical for the local re- 
modeling of axonal arbors and presynaptic sites regulating 
arbor size and shape in the vertebrate CNS. 

Materials and methods 

Zebrafish embryos 

All protocols were reviewed and approved by the Animal Care and Use Com- 
mittees of the RIKEN Brain Science Institute. WT embryos were obtained by 
natural matings, raised at 28.5°C, and staged by age and morphology 
(Kimmel et al., 1995). WT embryos were obtained from RIKEN WT, tupfel 
longfin (TL), and AB strains. TL's were provided by the Zebrafish International 
Resource Center (which is supported by grant P40 RRO 12546 from the Na- 
tional Institutes of Health-National Center for Research Resources) and by 
H. Baier and K. Slonchez of the Max Planck Institute for Neurobiology, 
Martinsried, Germany, asf'^''^^ (as^/robo2 mutant embryos were obtained 
by crosses of homozygous mutants, provided by C.-B. Chien and M.-Y. Low 
(University of Utah, Salt Lake City, UT), and asf'^^^'-/+ was generated by 
crossing to RIKEN WT or TL strains. 

In situ hybridization 

In situ hybridization was performed as described using digoxigenin-lobeled 
antisense probes (Schulte-Merker et al., 1 992; Lee et al., 2001 ) of casp-3a 
(a gift from M. Yamashita, National Research Institute of Fisheries Science, 
Yokohama, Japan; Yabu et al., 2001), casp-8 (GenBonk accession no. 
BC08 1 583; clone ID 235 1 947 obtained from Thermo Fisher Scientific; Eimon 
et al., 2006), and casp-9 (GenBank accession no. BC085434; obtained from 
the zebrafish gene collection, clone ID 7226798; Thermo Fisher Scientific) on 
80-hpf embryos. Embryos were fixed in 4% PFA for 2 h at RT or overnight at 
4°C followed by a graded methanol series to 1 00% methanol overnight at 
— 20°C. Embryos were rehydroted via a graded methanol series and washed 
3x for 5 min in PBST (PBS and Tween) followed by treatment with 0.1% 



collogenase for 50 min at RT. Embryos were then washed 3x for 5 min in PBST 
at RT followed by after fixation in 4% PFA for 20 min and washed for 3x for 
5 min in PBST before incubation for 60 min at RT with 2.5 (jl/ml acetic anhy- 
dride in 0.1 M triethanolomine. Embryos were washed for 2x for 10 min in 
PBST. Embryos were prehybridized for ~2 h at 65°C followed by hybridiza- 
tion with probe at 1 :50 or 1 : 1 00 dilutions overnight at 65°C. Incubation with 
10 [jg/ml RNoseA for 30 min at 65°C and two further 2x SSCT (SSC and 
Tween) washes. Reuse of sense and antisense probes helped to improve the 
signal to noise ratios. After probe hybridization and removal, embryos were 
washed 2x in 50% formamide and 2x SSCT for 30 min at 65°C followed by 
2x for 30 min each in 2x and 0.2x SSCT. Embryos were incubated 1-5 h 
with 2% blocking reagent (Roche) in maleic acid buffer before being incu- 
bated overnight at RT with a 1 :5,000 dilution of ontidigoxigenin alkaline 
phosphatase Fab fragments (Roche) in 2% blocking reagent. Unbound Fob 
fragments were removed with washes for 25 min in 2% blocking reagent and 
3x for 25 min in maleic acid buffer followed by incubation in BM purple 
(Roche) for colorimetric staining. The color reaction was stopped by 3 x 5-min 
PBST washes followed by fixation in 4%, PFA for 2 h at RT and transfer to 1 00% 
methanol as described earlier in this paragraph. After rehydration in PBST, 
clearing of the embryos was performed through a graded glycerol series to 
70% glycerol and stored at 4°C before imaging. Whole-mount in situ hybrid- 
ized embryos were embedded in agarose in PBS, fixed in 4% PFA, sectioned 
at 60 pm on a vibratome (DSK Microslicer DTK-3000W; Ted Pella), cover- 
slipped in 80% glycerol, imaged, and photographed using a microscope 
(Axioplan 2; Carl Zeiss) with a 40x Achroplan, NA 0.65 or 20x Achroplan, 
NA 0.5, a camera (AxioCom HRc; Carl Zeiss), and Axiovision (Carl Zeiss) or 
MetaMorph (Molecular Devices) software. 

DNA constructs and MO antisense oligonucleotides 

Constructs used in this study were as follows: isl2b {islet2b):GAL4VPl 6 
(originally named islet3:GAL4VPl 6; Campbell et al., 2007), the zebrafish 
islets gene was renamed isl2b (Pittmon et al., 2008), UAS:SCAT3 and UAS: 
SCAT3DEVG (control; gifts from M. Miura and Y. Yamaguchi, The University 
of Tokyo, Bunkyo, Tokyo, Japan; Tokemoto et al., 2003, 2007), and UAS: 
YFP-RabS (Campbell et al., 2007). The following MO antisense oligonucle- 
otides were obtained from Gene Tools: casp-3aM01, 5 '-TTGCGTCCACA- 
CACAGTCTCCGTTCAT-3', which has been shown to be effective in 
zebrafish (Yamashita et al., 2008); cosp-9M01, 5'-GAATCTGTCTGT- 
GTTTCTGCTCCAT-3'; casp-9SBM01 (Sidi etal., 2008), 5 '-GATGGAAAAA- 
CACACTTACGGACTG-3'; control (coMO; five-base pair mismatch for 
cosp-3aM01 ; lower coses letters represent the base changes), 5'-TTcCGTa- 
CACAgAGTCTCCcTTgAT-3'; s/Zt/atronsMOl , 5 '-GACAACATCCTCCTCTC- 
GCAGGCAT-3' (Barresi etal., 2005; Campbell etal., 2007); rofao2SBMO, 
5'-TAAAAAGTAGCGCAACTCACCATCC-3' (Wyott et al., 2010); casp- 
8M01 , 5'-ACAGGGTTTTAACTCACAGTAGATC-3' (Sidi et al., 2008); and 
p38mapkMO], 5'-AGTGGGTCTTTCTTTCTGCGACATG-3' (Hsu etal., 
201 1). robo2SBMO and s//t/a transMOl were injected at 2- or 1-ng/nl 
doses, respectively, which did not affect overall arbor growth on their own, 
except for the slitTa transMOl experiments described in Fig. 5 (E-H) in 
which 2 ng MO was injected, which increases arbor size and complexity 
(Campbell et al., 2007). For the epistosis experiments examining arbor dy- 
namics in Fig. 8, combinations of doses of MOs were used that on their own 
did not affect arbor dynamics (s//7 Jo transMOl [low dose, 1 ng], rotio2SBMO 
[low dose, 2 ng], and casp-3aM01 [low dose, 0.5 ng]). MOs were diluted 
(combined with DNA constructs) and injected to a final volume of ~1 nl, 
being one fifth the size of the yolk, into one- to four-cell stage embryos. MOs 
against the target genes ore effective in zebrafish, and because high doses 
of casp-3oM01 , casp-8M01 , casp-9M01 , cosp-9SBM01 , s/Zf/o transMOl , 
rofao2SBMO, and p38mapkMO\ con generate morphological phenotypes, 
the doses were chosen to avoid these and represent a partial loss of function 
(Campbell et al., 2007; Sidi et al., 2008; Yamashita et al., 2008; Wyott 
etal., 2010; Hsu etal., 2011). 

Labeling of RGCs and tectol cells 

Single RGCs or tectal cells were labeled by mosaic expression after coinjec- 
tion into one to four cell embryos of the plosmids isl2b:GAL4VP 16 and UAS: 
SCATS, UAS:SCAT3DEVG (control), or UAS:YFP-Rab3. Embryos were 
screened at ~76 hpf for one or a few RGC axonal arbors in the coudome- 
diol region of the optic tectum. This area is relatively flat, in the center of the 
tectum and where RGC arbors ore close to the surface of the brain, avoiding 
the likelihood of being obscured by the autofluorescence of the skin (Hutson 
et al., 2004; Campbell et al., 2007). Arbors that overlapped with or were in 
close proximity to other arbors, making it difficult to distinguish individual ar- 
bors, were eliminated from our analyses. RGC cell bodies for imaging were 
located in the central retina. Tectal cell bodies for imaging were located in 



BBB JCB • VOLUME SD3 • NUMBER 4 • 201 3 



the periventricular region of the optic tectum, close to the surface of the brain. 
Because SCATS or control construct expression level varied and to minimize 
any effects as a result of arbor morphology or movement on the FRET ratio 
between RGC arbors, we examined a brood population (19-28) of arbors 
per condition and observed consistent differences in the Venus/ECFP ratio 
between SCATS and control constructs. The effect of expression levels on the 
Venus/ECFP ratio appeared minimal. WT- and SCATS-expressing RGC ar- 
bors were indistinguishable in arbor size and branch tip dynamics from those 
expressing YFP-RabS. Arbor parameters were as follows: branch tip number, 
WT YFP-RabS 4.2 ± 0.67 compared with SCATS 4.086 ± 0.39 (P > 0.1); 
arbor area, WT YFP-RabS 150.60 ± S9.70 pm^ compared with SCATS 
1 57.7S ± 7.45 pm^ (P > 0. 1 ); and total arbor length, WT YFP-RabS 46.64 ± 
1 0.08 (jm compared with SCATS 45.26 ± 1 .90 (jm (P > 0. 1 ). Arbor dynam- 
ics were as follows: branch tip additions, WT YFP-RabS 737.19 ± 65.2% 
compared with SCATS 793.5 ± 96.14% (P > 0.1); branch tip retractions, 
WT YFP-RabS 752.57 ±131% compared with SCATS 778.74 ±131 .79% 
(P > 0.1); and mean total branch tip lifetimes, WT YFP-RabS 55.86 ± 
5.3 min compared with SCATS 54.95 ± 4.01% (P > 0.1). 

SCAT3 reporter imaging of RGCs and tectal cells and image analysis 

Imaging of RGCs or tectal cells expressing SCATS or control constructs was 
performed using a confocal microscope (FVIOOOD [Olympus] or LSM 780 
[Carl Zeiss]). Embryos were restrained by onoesthetization in 0.004% tric- 
aine and mounted dorsal side down in 1 .5% low melting point agarose 
(NuSieve; Lonza) in a Lab-Tek two-chambered #1 .0 borosilicate covergloss 
system (Thermo Fisher Scientific). Imaging was performed using 440- or 
458-nm laser excitation with 60x U Plan S Apochromat, NA 1 .2 water, 
60x U Plan S Apochromat, NA 1.3 silicone oil, or 63x C-Apochromot, 
NA 1 .2 water objectives. Emission was detected at 460-500 nm (ECFP) 
and 515-615 nm (Venus). 10-30 1- or 2-(jm z stacks were collected (512 x 
5 1 2 pixels). Care was token to avoid pixel saturation. For time-lapse imag- 
ing with SCATS, z stacks were collected every 5 min for ~2 h at ~26°C. 
The viability of the embryos was confirmed with a beating heart and blood 
flow after the imaging period. Only samples with minimum movement 
(where movement was not observed on images during confocal loser scan- 
ning) were analyzed, and individual time points were treated indepen- 
dently with their own internal reference of, for example, background 
fluorescence as described in the next paragraph. Stacks were aligned to 
eliminate x-y drift using the StackReg plugin in ImogeJ (National Institutes 
of Health). 

Images were analyzed as maximum intensity z projections of the 
ECFP and Venus channels using ImogeJ software. A stock of the Venus and 
ECFP z projections was created to precisely align the images. Venus/ECFP 
ratios over the whole RGC arbor or tectal cell body were determined using 
the polygon selection tool at o zoom of 200% to trace around oil the branch 
tips and intermodol regions of the arbor to the branch point closest to the 
parent axon or cell body on the Venus channel followed by calculating the 
mean fluorescence intensity per unit area for both the Venus and ECFP chan- 
nels. Mean background fluorescence intensity from identical regions of each 
channel was subtracted from the mean fluorescence intensity measurements, 
the ratio calculated and expressed as o percentage of Venus/ECFP Back- 
ground-subtracted Venus/ECFP ratios at branch points and surrounding re- 
gions were calculated as earlier in this paragraph except using the polygon 
selection tool to trace around branch points and three surrounding regions 
on the same arbor, adjacent to the branch point as described in Fig. 1 (B-D 
and F-H). The mean Venus/ECFP ratio of the three surrounding regions is 
depicted in the graphs (Figs. 1 T, 2 F, and 5 Gj. Representative intensity 
modulated display (IMD) images were generated using the ratio images 
function in MetaMorph software. Measurement of arbor parameters on the 
YFP and IMD ratio images was very similar, showing that the whole of the 
arbor is visible in the ratio images (Fig. S2, A-D). Images presented in 
Figs. 1 , 2, 3, 5, 6, 7, S3, S4, and S5 and Videos 1 , 2, 3, 4, 5, 6, 7, 8, 9, 
and 10 were median filtered (radius of 1 pixel) to reduce noise, and input 
levels were set to 1 00 or 1 50 to increase brightness for displaying. Time in- 
dicators were added to the Videos 1 , 2, 3, 4, 5, 6, 7, 8, 9, and 10 using 
the Time Stamper function, and arrowheads in Video 1 were added with the 
Arrow Maker Tool plugin in ImogeJ. 

Branch tip and presynaptic puncta time-lapse imaging 
and image analysis 

Imaging of arbors expressing YFP-RabS was performed using on upright 
(LSM 510 META; Carl Zeiss), inverted (LSM 780), or an upright confocal 
microscope (FVIOOOD) as described previously (Campbell et ol., 2007). 
YFP-RabS was excited with 515- or 514-nm laser lines, and emission was 
detected with o 530-600-nm band pass filter or o BA535-564-nm barrier 



filter. 10-35 2.0-pm z stacks (512 x 512 pixels) were collected every 
15 min for 5 h, incubated at 26°C with a 6Sx N-Achroplon, NA 0.9 long 
working distance water immersion objective on the LSM 510 META, a 6Sx 
C-Apochromot, NA 1 .2 objective on the LSM 780, or a 60x LLJMPFL, NA 
0.9 long working distance water immersion objective on the FVIOOOD. 
Care was token to ovoid pixel saturation. The viability of the embryos was 
verified after imaging via confirmation of a beating heart and blood flow. 
Time-lapse images were processed using the group z projector plugin for 
ImogeJ, and individual images were saved as an 8-bit TIFF stack. Because 
YFP-RabS is trafficked to both synaptic and nonsynoptic membrane com- 
partments, it is possible to use it to view both presynaptic sites and oxonal 
arbor morphology. Branch tip tracking was performed on unthresholded im- 
ages. The numbers and lifetimes of YFP-RabS-expressing presynaptic sites 
were quantified using thresholding in ImageJ software. A YFP-RabS punctum 
was defined as an area of >2 pixels whose intensity was >5x that of a non- 
punctate region of the some arbor and saved as an 8-bit TIFF stock 
(Campbell et ol., 2007). The fluorescence intensity of each arbor was nor- 
malized to a nonpunctate region of the some arbor, i.e., the parent axon. 
Tracking of branch tips and presynaptic puncta through different time points 
was performed manually using the ObjectJ for ImageJ plugin written by 
N. Vischer (University of Amsterdam, Amsterdam, Netherlands; source code 
con be found in Supplemental material). Total branch tip and presynaptic 
puncta lifetimes were defined as a branch tip or presynaptic punctum being 
present in one or more time-lapse imaging frames, calculated by summing 
up the number of frames that either was present for, and multiplied by 1 5 (the 
duration between frames) to generate the lifetime in minutes. 

Quantification of arbor parameters 

Maximum intensity projections of YFP-RabS-expressing arbors were quan- 
tified OS in Campbell et al. (2007): The number of branch tips >4 pm are 
counted. The total arbor length is the sum of all branch tip lengths and in- 
ternodol distances. The arbor area is defined by a convex polygon joining 
all branch tips. The order of an arbor is the mean number of branch points 
counted between a branch tip and the first branch point of the arbor. For 
puncta number, images were thresholded as in the preceding paragraph 
and counted using the analyze particle function with a pixel size of >2 in 
ImogeJ software. 

Caspase-3 activation of SCATS 

Tectal cell bodies of 76-80-hpf zebrafish larvae expressing SCAT3 or control 
construct (control) in one or a few tectal cells were imaged on an inverted 
FVIOOOD confocal microscope as described in the section SCATS reporter 
imaging of RGCs and tectal cells and image analysis. After imaging, embryos 
were remounted (dorsal up) in 1 .5% low melting point agarose (NuSieve) 
on a gloss microscope slide. The Cospose-S activator PAC-1 (procaspose- 
octivating compound 1 ; Tocris Bioscience; Putt et ol., 2006) was dissolved in 
DMSO followed by embryo medium to a final concentration of 40 pm and 
pressure injected into the tectal ventricle, an approach previously shown to be 
able to target reagents to tectal cells (Leu and Schmidt, 2008). Embryos were 
placed individually in 500 pi embryo medium in a 24-well plate and incu- 
bated at RT for 3 h, sufficient time for PAC-1 to induce a high proportion of 
TLJNEL-positive apoptotic cells in the retina and tectum (described in the next 
section). Embryos were subsequently remounted for an inverted confocal 
microscope (FVIOOOD), and imaging and analysis were performed as de- 
scribed earlier in this paragraph. 

TUNEL labeling 

TUNEL labeling was performed using the In Situ Cell Death Detection kit, 
polymerized horse radish peroxidase (Roche). 80-hpf embryos were fixed in 
4% PFA and permeobilized with 1 00% methanol at -20°C followed by 0. 1 % 
collagenose at RT and PBST washes. Samples were incubated in TUNEL reac- 
tion mixture for 60 min at 37°C in the dork. After PBST washes, samples were 
incubated with antifluorescein-polymerized horse radish peroxidase (Roche) 
at a dilution of 1:1,000 overnight, washed in PBST, and stained with 
0.5 mg/ml S,3-diaminobenzidine with 0.0003% H2O2 in PBST 

Embryos were mounted in 70% glycerol and placed in a gloss plate 
and imaged using a stereo dissecting microscope (SZX16; Olympus) with 
o super depth of focus Plan Apochromat 1 x parfocal objective zoom 5 and 
0 camera (DP50-CU; Olympus) at 1,392 x 1,040 pixels with Studio Lite 
software (Olympus). Preincubation with 20 pm of the Bcl-2 (B cell lym- 
phoma 2) inhibitor HA14-1 (Sigma-Aldrich) for 1 h was used as o positive 
control for TUNEL labeling. 

Acceptor photobleaching of SCAT3 

To test for the presence of FRET between the ECFP and Venus moieties of 
SCATS and the uncleovoble control construct (control), 80-hpf RGC arbors 
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and tectal cell bodies were imaged before and after partial photobleach- 
ing by exposure to 1 00% laser power of a 5 1 4-nm loser line from a multi- 
argon line laser (Melles Griot) for 5 min. Reimaging of the Venus/ECFP 
ratios of ttie same cells was performed as described in the section SCAT3 
reporter imaging of RGCs and tectal cells and image analysis. 

Imaging RGC pathfinding to the optic tectum 

For lateral views, embryos with YFP-Rab3-labeled RGC axons were processed 
for anti-GFP immunohistochemistry. Approximately 80-hpf embryos were fixed 
and permeobiiized as described for in situ hybridization. After the collagenase 
PBST washes, embryos were blocked in 1 0% fetal calf serum for 1 h at RT fol- 
lowed by incubation overnight at 4°C with a chicken primary anti-GFP poly- 
clonal antibody (Aves Labs, Inc.) at a dilution of 1 :3,000. The primary 
antibody was washed off with 4 x 30-min washes in PBST at RT followed by 
incubation for 4 h with on anti-chicken Alexo Fluor 488-coniugated sec- 
ondary antibody (Molecular Probes and Invitrogen), diluted 1 :500. After 4 x 
30-min PBST washes, embryos were mounted in 1 .5% low melting point aga- 
rose for an upright confocol microscope (510 META). Images were captured 
as a z series (step sizes of 2 pm) of 5 1 3 x 5 1 2-pixel scons using a long work- 
ing distance 20x, NA 0.5 N-Achroplan water immersion objective and a digi- 
tal zoom of 0.8. 

Statistical analyses 

Statistical analyses were performed using Prism version 4 software (GrophPod 
Software). As it con be difficult to determine accurately whether dotosets fol- 
low a normal Gaussian distribution or not, we applied nonporometric statisti- 
cal analyses, which ore usually more stringent, using the Monn-Whitney 
U test, allowing for whether standard deviations of the conditions being com- 
pared were significantly different when comparing two conditions. When 
comparing three or more conditions, the Kruskal-Walloce test with a post hoc 
Dunn's multiple comparison test was used. P < 0.05 was considered statisti- 
cally significant. All results ore presented as mean values ± SEM. 

Online supplemental material 

Fig. SI shows cospases ore expressed in RGCs and the tectum. Fig. S2 shows 
that SCATS reports Caspase-3 activation. Fig. S3 shows cospase activation 
is undetectable in tectal dendrites. Fig. S4 shows caspase activation remains 
lower after branch tip retraction. Fig. S5 shows that in caspase morphants, 
RGC axons follow normal paths to the optic tectum. Video 1 shows the 
dynamics of caspase activity in young RGC arbors. Video 2 shows the dy- 
namics of caspase activation in older arbors. Video 3 shows that WT arbors 
ore highly dynamic. Video 4 shows that coMO arbors are highly dynamic. 
Video 5 shows that casp-3a knockdown leads to more stable arbors. Video 6 
shows that casp-9 knockdown leads to more stable arbors. Video 7 shows 
that p38mapk knockdown leads to more stable arbors. Video 8 shows 
that low dose MO against casp-3a does not affect arbor dynamics. Video 9 
shows that low dose MO against slitla does not affect arbor dynamics. 
Video 10 shows that a combination of low dose casp-3oM01 and low 
dose slitla tronsMOl leads to more stable arbors. A source code for 
Object] plugin for ImogeJ is also provided in a ZIP file (containing source 
code ob|ectj_-o.96d.jar, branch macro branches-2.txt, branches object file 
Branches-2.oii, and a demo file demo-1 .tiff). Online supplemental material is 
available at http://www.jcb.org/cgi/content/full/jcb.201 303072/DCl . 
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